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Abstract. In this paper, we propose a parameterization which is based on a phenomenological model 
involving the wounded quarks interactions for explaining the average charged particle multiplicity 
{n c h) and the central pseudo-rapidity density {(dn/dri) rl ^o) and their variations with the center-of-mass 
energies (y^ijvjv) in relativistic heavy-ion collider experiments. The model also interrelates nucleus-nucleus 
(A-A) collisions with p-A and p-p interactions. Our parameterization rests on simple assumptions 
regarding mean number of participating quarks and their average number of collisions. The results 
for {n c h) and their variations with the mass number of colliding nuclei and ^snn are well supported 
by the experimental data. The results obtained from modified Glauber model involving wounded 
nucleons are also shown for comparison. Finally, we give the predictions from our model for A-A 
collisions at the Large Hadron Collider (LHC) and Compressed Baryonic Matter (CBM) experiments. Our 
results indicate towards a possible universal production mechanism existent for p-p, p-A and A-A collisions. 
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1 Introduction 

The ultimate goal of ultra-relativistic heavy-ion collision 
experiments is to test the predictions of quantum chro- 
modynamics (QCD) because these can reveal the nature 
of hadronic interactions at extreme temperature and/or 
density and throw light on the role played by the quarks 
and gluons in multiparticle production mechanism [1] . In 
order to verify the existence of a deconfmcd plasma cre- 
ated in nuclear collisions, we require accurate and precise 
knowledge about its subsequent hadronization. Hadronic 
multiplicities and their correlations can reveal enough in- 
formation on the nature, composition and size of the fire- 
ball formed in the collider experiments [UGH ■ The evolution 
of the measured global observables in such experiments 
varies with the size and mass of the colliding nuclei and it 
emphasizes that significant rescatterings or multiple scat- 
terings occur and contribute to the evolution of the energy 
deposition in the fireball. 

The search for some regularities and systematics present 
in the average multiplicity distributions of hadron-hadron 
(h-h), hadron- nucleus (h-A) and nucleus-nucleus (A-B) 
collisions has grown into a fascinating topic in strong inter- 
action physics today because these properties hint at the 
underlying universal production mechanism for charged 
particles. The study of p — p and p — A collisions is par- 
ticularly useful because it helps our understanding, e.g., 
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in which way A — B collisions differ from simple super- 
position of nucleon-nucleon collisions. However, our un- 
derstanding of the subject is still quite poor partially be- 
cause of the reason that in soft hadronic or nuclear col- 
lisions, the role of perturbative QCD is less clear as it 
fails close to the phase hadronization point where non- 
perturbative effects become clearly dominant. Several new 
experimental information on the multiparticle production 
have been reported in the recent past. Consequently a lot 
of efforts have also been put forward to understand and/or 
organise the experimental data by using various theoret- 
ical as well as phenomenological schemes [3,5,6,7,8,£)l[TUl 
P1P^P^pll[T5l[T6l[r7lP^p^[^[^[^[^[Ml[^[^ . Some 
of these models explain multiplicities in nucleus-nucleus 
collisions mostly by using either quark-quark interactions 
or by considering basic nucleon-nucleon scatterings [23, 
[MMM2H]- However, predictions of models that success- 
fully describe multiparticle production at RHIC energy 
usually differ at LHC energy by a factor of two [27] which 
signifies a major discrepency. So we focus our efforts in 
formulating a model which can describe multiplicity dis- 
tributions consistently in the entire energy range accesible 
today at Relativistic Heavy Ion Collider (RHIC), Large 
Hadron Collider (LHC) and/or in future at Compressed 
Baryonic Matter (CBM) experiment of the FAIR (Facility 
of Antiproton and Ion Research) at GSI Darmstadt. 

How one can describe the mechanics of such a relativis- 
tic many body problem in terms of simple, basic physical 
principles pose a major challenge in high energy heavy- 
ion physics today. Nucleus-nucleus collisions differ from 
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nucleon-nucleon collisions in many respects. The presence 
of many constituent quarks participating in a chain of col- 
lisions makes the picture too complicated. In this paper, 
we describe nucleus-nucleus collisions in terms of effective 
number of constituent quarks participating in the colli- 
sions and the effective number of collisions suffered by 
each of them. In this model, we propose a gluon being ex- 
changed between a quark of projectile or first nucleus and 
a quark belonging to target or other colliding nucleus. The 
resulting colour force is thus somewhat stretched between 
them and other constituent quarks because they try to re- 
store the colour singlet behaviour. When two quarks sep- 
arate, the colour force builds up a held between them and 
as the energy in the colour held increases, the colour tubes 
thus formed finally break-up into new hadrons and/or 
quark- antiquark pairs. We have essentially modihed the 
phenomenological model earlier proposed by Singh et al 
[25] . We consider a multiple collision scheme in which a 
valence quark of the incident nucleon suffers one or more 
inelastic collisions with the quarks of target nucleons. In 
this process, quark loses energy and momentum which re- 
sults into creation of new hadrons. Further, we assume 
that all the quark collisions are independent and their ef- 
fects are incoherently superimposed. We hnd that charged 
particle multiplicities arising from p — p, p — A and A — B 
collisions can be accounted in a consistent manner by ade- 
quately incorporating the number of participating quarks 
and mean number of quark-quark collisions. The rest of 
the paper is organised as follows. In section II, the de- 
tailed description of the model is presented, the results 
and discussions are summarized in section III and finally 
in section IV, main conclusions and future prospects of 
this study are succinctly given. 



2 Description of Model 

If we assume a universal mechanism of charged parti- 
cle production in the hadron-hadron, hadron-nucleus and 
nucleus-nucleus collisions, it must be driven by the avail- 
able amount of energy involved for the secondary produc- 
tion, and it must depend on the mean number of partic- 
ipant quarks |2S|. The main ingredients of our model are 
taken from the paper by Singh et. al. [28 . There exist 
several papers in the literature emphasizing various types 
of fitting parametrizations in order to provide a unified 
description for the produced charged particles in p — p col- 
lisions at various energies. We have taken some of these 
functions with the values of their constants from the liter- 
ature 30 , 31 , 32 . 33 and fitted them with the experimental 
data [3^[331[3^||371|551[3^1|iTJl|iT||?2] as shown in Fig. 1. All 
these parametrizations fit the experimental data for p — p 
collisions at intermediate and higher energies but usually 
show some disagreement with the experimantal data at 
lower energies ( see Fig. 1). Recently, Jeon and collabora- 
tors have shown that the total mean multiplicity obtained 
at RHIC can be bounded by a cubic logarithmic term 
in the energy [43] . Based on the findings of Ref. [43] , we 
propose here a new parameterization which extends upto 
a cubic logarithmic term to accomodate the entire p — p 



experimental data [MUnSllMlEIlEHJIMllinillllllS starting 
from low energies (i.e. from 6.15 GeV) upto the highest 
LHC energy (i.e. 7 TeV) as follows : 



< n ch > hp = (a'+b'lny/s^+c'ln 2 ^/s^+d'ln 3 ^/s^)- 



(1) 

In Eq. (1), a is the leading particle effect and ^fs2 is the 
available center-of-mass energy (i.e., = ^fs — tub — 
m T [33], where is the mass of projectile and tut the 
mass of the target nucleon, respectively), a', &', c' and d! 



are chosen constants. We take the values 



1.8, b' 



0.37, d = 0.43 and d' = 0.04 as derived from the best 
fit to the data. The value of a is taken as 0.85 in our 
calculations [25] . 

We can extrapolate the validity of this parametrization 
for the produced charged particles in the hadron-nucleus 
interaction by considering the multiple collision effect in 
the nucleus. The number of constituent quarks which par- 
ticipated in hadron-nucleus (h-A) collisions share the total 
available center-of-mass energy ^Jsua and thus the energy 
available for each interacting quark becomes <Jsh~A~/N q %A [28] 
where N q A is the mean number of constituent quarks 
in h-A collisions. The total available squared center-of- 
mass energy ShA in h-A collisions can be related to s a as 
ShA — v q A s a [34,45 with v q A as the mean number of 
inelastic collisions of quarks with target nucleus of atomic 
mass A. Within the framework of the Additive Quark 
Model 46 , 47 , 48 ,49,50], the mean number of collisions in 
hadron-nucleus interactions is discussed in relation to the 
ability of constituent quarks in the projectile hadron to in- 
teract repeatedly inside a nucleus. The average number of 
constituent quarks (N^ A ) interacting inelastically in h-A 
collisions is given by |51j . 



VqA ' 



(2) 



where, v(A) is the mean number of collisions of the inci- 
dent hadron and is equivalent to [STj . 



u(A) = v qA - 



qA 



(3) 



'hA 



with a™ A as the inelastic hadron-nucleus cross-section and 
N c is the number of valence quarks, D q A represents the 
mean number of collisions of the wounded quark inside 
the nucleus A and is defined by [STj . 



An 



VqA 



qN 



(4) 



where cr q n N and a l qA are the inelastic cross-sections for 
quark-nucleon (q-N) and quark-nucleus (q-A) interactions, 
respectively and A is the atomic mass of the target nu- 
cleus. 

Now, using the Eq. (2), (3) and (4), the average num- 
ber of constituent quarks (N q ) in h-A collisions can be 
rewritten as [Z51T51"] . 



N hA = N ^ qA 
l\ q 



a 



(5) 



hA 
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The quark-nucleus inelastic cross-section (cr*^) is deter- 
mined from o™ N (w by using Glauber's approxi- 
mation given as follows |28j : 



AB 



7 qA 



d 2 b 



(6) 



where profile function D A {b) is related to the nuclear den- 
sity, p(b, z) by the relation [35] : 



D A (b) 



p(b, z)dz, 



(7) 



and we have chosen the form for the nuclear density as 
follows [S2]: 



P(b, 



Po 



1 + exp{ 



(8) 



where all the notations have their usual meaning 1521 . We 
have taken the values of the constants i.e., po, R and a for 
different colliding nuclei separately from the Ref. [53] . 

By using Eq. (5), one can easily see that the average 
number of constituent quark (N q ) is 1 for h-p collisions 
as in nucleon-nucleon collisions the interactions occur be- 
tween single constituent or dressed quark in accordance 
with the additive quark picture [4"9"115T)] and other quarks 
are considered to be spectators. Thus the energy of the 
initial thermalized state which is responsible for the pro- 
duction of secondary particles is mainly due to the inter- 
acting single quark pair. The spectator quarks which are 
not the part of thermalized volume at the moment of col- 
lision, do not participate in secondary particle production 
process. As a result, the leading particles [54"1I5"5] resulting 
from the spectator quarks carry away a significant part of 
the energy [5^ . Finally, the expression for average charged 
hadron multiplicity in h-A collisions can be expressed as 
follows [551: 



< rich >hA= N. 



hA 



MhA 

q 



+ c'ln 2 



MhA 
Q 



d'ln 3 



MhA 

q 



(9) 



The generalization of the above picture for the case of 
nucleus- nucleus collisions can be achieved as follows 1251: 



< rich >ab= N, 



AB 



a' + b'ln 



y/SAB 

N AB 
q 



Furthermore, mean number of participating quarks N t 
in A-B collisions can be calculated by generalizing Eq. (8) 
in the following manner |28) : 



N AB = t 

q 2 



>qA 



N A a 



qB 



'AB 



'AB 



(12) 



where is the inelastic cross-section for nucleus-nucleus 
(A-B) collision and can be expressed in the following man- 
ner [2SJI3B]: 



a AB 



AV3 + B l/3 



,41/3 + £1/3 



(13) 



where c is a constant and has a value 4.45 for nucleus- 
nucleus collisions. One can see that the extended parametriza- 
tion in Eq. (10) relates nucleus-nucleus collisions to hadron- 
nucleus and hadron-proton collisions and the values of the 
parameters a', b' , c', and d' remain unaltered which shows 
a universality in the role of basic quark-gluon interactions 
in all these processess together. 

In a search for creating quark gluon plasma (QGP), 
greater emphasis is laid on the central or head-on col- 
lisions of two nuclei. The resulting mean multiplicity in 
central collisions can be obtained by using Eq. (10) in the 
following manner : 



< n ch 



•^cejitral _ 



+ b'ln(K B ^) 1/2 +c'ln 2 (^ B s a )^ 2 



+ d>ln 3 (^B Sa) i/2 



(14) 



(10) 



where ^Js A b = A \j v q B s a [28 and the mean number 
of inelastic quark collisions v^ B is [28] : 



The pseudorapidity distribution of charged particles is 
another very important quantity in the studies of particle 
production mechanism in high energy h — h and A — B 
collisions, which, however, is still not understood properly. 
It has been pointed out that (dn c h/dr]) can be used to get 
the information on the temperature (T) and density (p) 
of the QGP [5711551159] . Moreover, mid-rapidity density 
reflects the different stages of the reactions. To calculate 
the pseudorapidity density of charged hadrons, we first fit 
the experimental data of (dn c h / 'drj)^fL Q by using a sim- 
ilar parametrization as used earlier in Eq. (1), but with 
slightly different values of constants. Further, we calculate 
the central pseudo-rapidity density (dn c h / drj)^£ Q by using 
the basic idea of two component model for nuclear coll- 
sions [BD] , which suggests that the pseudo-rapidity density 
for A—B collisions comes from two parts i.e., contributions 
arising due to hard processes and the other due to the soft 
processes. The contributions from hard processes is propo- 
tional to the total number of collisions i.e., N^ B u^ B and 
the portion from soft processes is propotional to the par- 
ticipating quark constituents i.e., N^ B only Thus, the ex- 
pression for (dn c h/dri)^£ Q can be parameterized in terms 
of p-p rapidity density as follows : 



AB 



VqAVqB 



A(J qN B(J qN 

-in ' -in 
a qA G qB 



(11) 



dn, 



ch 



dr) 



AB 



t;=0 



/ dn ch 
\ dr] 



r?=0 



[(l-x)Nj 



AB 



xNt B vt B \ 



(15) 



4 



Ashwini Kumar et al.: Wounded Quarks and Multiplicity at Relativistic Ion Colliders 



■ INELdata 

□ NSD data 

Our parameterization Eq. (1) 

a, + b,lns + Cjhrs 

a, + b,s 1/4 

a, + b,s" 




Fig. 1. Variations of total mean multiplicities of charged 
hadrons in p — p collisions with ^/snn- 
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Fig. 3. Variations of total mean multiplicities of charged 
hadrons in central collisions for different colliding nuclei with 



Modified Glauber model 
Our Model Results 
Experimental data 



^(GeV) 

Fig. 2. Variations of total mean multiplicities of charged 
hadrons in central Au — Au collisions with ^Jsnn- Experimen- 
tal data points are taken from Ref. 62 , 63 : 64 65 . 



where x quantifies the relative contributions of hard and 
soft processes. The fraction x corresponds to the hard pro- 
cesses and the remaining fraction (1 — x) thus describes 
the contributions from the soft processes. In Eq. (15), we 
have used the parameterization for (dn c h / drff^L^ same as 



used in Eq. (1) for calculating < 
and can be given as follows: 



rL c h > in p-p collisions 



< {dn c h/dr))^ > — (a 1 + b 1 ln y /s£+ c^n 2 

+ d'^n 3 ^/!^) - a, (16) 

where a[ = 1.78, b[ = 0.025, c[ = 0.041 and d[ = 0.0017 
as derived from the best fit to the data. We take the value 
of x at different yfs NN according to Ref. [50] • 

3 Results and Discussions 

In Fig. (1), we present the inelastic (filled symbols) and 
non-single diffractive (NSD) data (open symbols) of charged 



hadron multiplicity in full phase space for p — p collisions 
at various center-of-mass energies from different experi- 
ments e.g., ISR, UA5 and E735 ^M^H^TMS^^SM 
|4"2"] . We use inelastic data at very low energies (filled sym- 
bols) because NSD data are not available for these ener- 
gies and also the trend shows that the difference between 
inelastc and NSD data is very small at lower energies. 
Further, we fit this data set with four different functional 
forms. The short-dashed line has the functional form as : 
a 2 + b 2 s 1 ^ 4 which is actually inspired by the Fermi-Landau 
model [3"0"ll3"Tll3"2"] . It provides a resonable fit to the data at 
higher y/sJ^N with a 2 = 5.774 and b 2 = 0.948 How- 
ever, since a 2 summarizes the leading particle effect, its 
value should not lie much larger than two. The dotted line 
represents the functional form as : a\+b± lns + ci ln 2 s |33] 
and it fits the data well at higher ^/snn but shows a large 
disagreement with the experimental data at lower center- 
of-mass energies. The dashed-dotted line represents the 
form a 3 + 6 3 s™ [33] and it also provides a good qualita- 
tive description of the data with a 3 = 0, 63 = 3.102 and 
n = 0.178 [40] . The solid line represents our parametriza- 
tion given by Eq. (1) and it undoubtedly offers the most 
reasonable fit to the data starting from very low energy 
upto very high ^snn- 

Fig. 2, shows the variation of mean multiplicity of 
charged hadrons produced in central Au — Au collision 
with respect to ^/saT/v- We also compare our model results 
(solid line) and the modified Glauber model results (dash- 
dotted line) [61] with the experimental data of AGS and 
RHIC 62,63,64,65]. We would like to mention here the 
main difference between modified Glauber model and the 
standard Glauber calculation. In standard Glauber model, 
a participating nucleon can have one or more number of 
collision at the same v^tvat an d with the same value of 



(V* 



NN/ 



However, in modified Glauber model, a par- 
ticipating nucleon loses a fraction of its momentum after 
the first collision and participate in subsequent collisions 
with somewhat lower energy and also at a lower value of 
a inei(V~s)- This modification significantly suppresses the 
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Table 1. Total mean multiplicity of charged hadrons produced 
in d — Au collision at ^Jsnn = 200 GeV. 





Colliding Nuclei (n ch ) mm - btas 


/ \ central 

{rich) 


Our Model 
PHOBOS E2 


d- Au 82 
d - Au 87ig 


162 




Fig. 4. Variations of pseudorapidity density at mid-rapidity 
for p — p collisions with ^/snn- Filled symbols are data from 
inelastic p — p events and open symbols are experimental data 
for NSD events [35H36U37U38II42II68II69II701I7TU72] . Solid line is 
the outcome of our parametrization. 



number of collisions (N co u ) in comparison to that obtained 
in the standard Glauber model. We find that the results 
obtained in our model give a good description to the ex- 
perimental data in comparison to the modified Glauber 
model predictions, especially at lower center-of-mass en- 
ergies. 

In Fig. 3, we show the variation of mean multiplic- 
ity calculated from our model for Au — Au, Pb — Pb and 
Cu — Cu central collisions with respect to ^snn- We also 
compare our model results with the experimantal data of 
RHIC and SPS experiments [62,63,64,65,66 . One can see 
from Fig. 3, that the model results give excellent fit to the 
separate experimental data points for Cu — Cu, Au — Au 
and Pb — Pb, respectively. This concludes that our phe- 
nomenological model works well because it considers the 
picture of incoherent superposition of basic quark-quark 
collisions in describing the nucleus-nucleus collisions. 

In Table I, we have shown the mean multiplicity of 
charged hadrons in minimum-bias as well as in central 
d — Au collisions at ^/sjvjv = 200 GeV as obtained in our 
model. We have also shown the experimental results 67 
for this collision for comparison. The good agreement be- 
tween the experimental and model results shows that it 
can be used for asymmetric colliding nuclei as well. 

In Fig. 4, we present the inelastic (filled symbols) and 
non-single diffractive (NSD) data (open symbols) of 
at mid-rapidity for p — p collisions at various center of 
mass energies from different experiments e.g., ISR, UA5, 
E735, RHIC and LHC ^^^MMMWiilMinM- 
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Fig. 5. Variations of pseudorapidity density at mid-rapidity 
for A — A collision with y / sjvjv. Symbols are data from various 
heavy ion collision experiments [62 , 63 64 , 65 66 , 73ll74| . 

Table 2. Predictions of total mean multiplicity of charged 
hadrons produced in heavy-ion collision experiments. 



Experiment 



LHC 



CBM I 



CBM II 



Colliding Nuclei 

^/snn (in GeV) 

(KchYourZldel 

(n ck ) ce * tra ' ibttsi 
K)ggS"' EE! 



Pb-Pb 
2760 
8895 
14377 
13124 (N part = 386) 



Au — Au Au - Au 

3.97 (E u , b ss 8.4) 8.1 {E, ab m 35A) 
328.75 670.8 
380.8 837 



; 308 



: 650 



Further, we fit (solid line) this data set with our functional 
form of Eq. (16), which is similar to Eq. (1) and obtained 
the values of constants as a x — 1.78, b 1 — 0.025, c x = 0.041 
and d 1 = 0.0017. The central rapidity hadron density thus 
continues to increase with the colliding energy at LHC by 
an added triple logarithmic form. 

In Fig. 5, we plot the variations of {dn c h/ drfj^Q ob- 
tained in our model for Au — Au, Pb — Pb and Cu — Cu 
central collisions with respect to We also compare 

our model results with the experimental data of RHIC, 
SPS and LHC experiments [^IIMIIMIISSIIMIIZIIZI] ■ Our 
model results compare well with the experimental data 
from very low energy to very high energy. As we didn't 
take any effect of final-state interactions in our model, it 
hints at the fact that the hadron multiplicity in A—B colli- 
sions is mainly driven by initial parton production and the 
effect of final-state interaction should be negligibly small 
as suggested in some earlier calculations [751F76] . 

In Table II, we give the predictions for the mean mul- 
tiplicity of charged hadrons, produced in minimum-bias 
events as well as in central events, based on our model 
calculations at LHC and CBM energies. Specifically, we 
calculate the charged hadron multiplicity in Pb — Pb colli- 
sions at ^snn = 2.76 TeV for LHC and the mean multi- 
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plicity of charged hadrons in Au—Au collisions at E^i = 8 
AGeV and 35 AGeV for future CBM experiment. We also 
compare our model results with the other model predic- 
tions like hadron string dynamics (HSD) model for CBM 
experiment [771175] and model of W. Busza for LHC ex- 
periment [6"1[79"]. respectively. 



4 Conclusions and Future Prospects 

Within the framework of constituent quarks, we have given 
a parametrization which interrelates the multiplicity and 
mid-rapidity density distributions in p-p, p-A and A-A 
collisions from a few GeV upto the highest LHC energies. 
It involves the interactions of the constituent quarks of 
the colliding objects or beam and target particle or nu- 
clei and the mean number of collisions suffered by them. 
We assume that the production of secondary particles is 
proportional to the fraction of the available energy for the 
participant quarks. This picture describes consistently the 
'soft' hadron production in p — p, p — p, p — A and A — A 
collisions. In heavy-ion collisions, one participating quark 
can have larger number of interactions due to the large 
size of the nucleus as well as due to a large travel path to 
be travelled inside the nucleus and hence, we get larger en- 
ergy available for the secondary particle production [25] , 
We find that our parametrization gives an excellent fit to 
the p — p data in the entire energy range. We also find 
that a suitable extension of this description works very 
well for the nucleus- nucleus (A — A) collisions in the entire 
energy range. We have also compared our model results 
with the results obtained from modified Glauber model 
to demonstrate the difference between the role played by 
nucleon-nucleon interactions vis-a-vis quark-quark inter- 
actions. We have also shown the charged particle mul- 
tiplicity in central as well as in minimum-bias d — Au 
collisions at ^/snn — 200 GeV and we find the results 
compare quite well with the experimental results. This ex- 
ercise clearly shows that our model also works quite well 
for asymmetric colliding nuclei case. We have also pro- 
posed an extension of the two-component model based on 
our parameterization and the results for rapidity distribu- 
tion obtained by us again work well. We have given our 
predictions for LHC and CBM data as well and compare 
them with other model predictions [51179) existing in the 
literature. 

In conclusion, we have attempted to point out a mech- 
anism for 'soft' hadron production common to hadron- 
hadron, hadron-nucleus and nucleus-nucleus collisions. Due 
to non perturbative nature, QCD is not expected to work 
in these cases. So lacking a workable theory, we have relied 
on the phenomenology and all the salient features of the 
experimental data are very well explained by our param- 
eterization. We hope that the work done here will high- 
light the multiparticle production mechanism present in 
nucleus-nucleus collisions in relation to p-p interactions. 
This can help in detecting any deviation observed in the 
data from the predictions of our model and thus a hint for 
QGP formation would appear likely in undergoing ultra 



relativistic heavy-ion experiments. Furthemore, future ex- 
tension of this scheme can help us to study the rapidity, 
pseudorapidity as well as transverse mass distributions of 
the charged hadrons in the entire rapidity region and this 
can form a matter for future publication. 
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